Nuclear spin driven quantum relaxation in LiYo.998Ho .oo2F4 
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Staircaselike hysteresis loops of the magnetization of a LiYo.998Hoo.oo2F4 single crystal are observed 
at subkelvin temperatures and low field sweep rates. This behavior results from quantum dynamics 
at avoided level crossings of the energy spectrum of sing le Ho 3+ ions in the presence of hyperfine 
interactions. Enhanced quantum relaxation in constant transverse fields allows the study of the 
relative magnitude of tunnel splittings. At faster sweep rates, nonequilibrated spin-phonon and 
spin-spin transitions, mediated by weak dipolar interactions, lead to magnetization oscillations and 
additional steps. 
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The problem of quantum dynamics of a two-level sys- 
tem coupled to an environment (boson or fermion bath) 
is at the core of mesoscopic physics We show that 
the new field of "mesoscopic magnetism" , which studies 
the tunneling of large magnetic moments in the presence 
of phonons and spins, is not limited to molecular com- 
plexes and nanoparticles, but it can be extended to other 
systems such as rare-earth ions. After the first studies on 
large spin molecules Mni2-ac |||| an d Fe§ § , the role of 
the spin bath on the tunnel mechanism was shown [^|-^) . 
In particular, quasistatic fields due to dipolar interactions 
between molecules lead to a distribution of internal fields, 
and field fluctuations, essentially of nuclear spins, give 
homogeneous level broadening allowing the restoration 
of tunneling in a finite energy window, at low tempera- 
ture; this broadening being much larger than the phonon 
one, it is more relevant to induce tunneling. This mech- 
anism is efficient unless all nuclear spins of the molecule 
are frozen, which occurs only below the mK scale. In low 
spin molecules, large tunneling gaps favor spin-phonon 
transitions. Although the hyperfine induced level broad- 
ening is the same as in large spin molecules, the phonon 
bath becomes as important as the spin bath p0|| . In all 
these cases, the role of field fluctuations was clearly evi- 
denced. 

This description is for the relatively weak hyperfine 
interactions of Mni2 or Fes molecules, and therefore for 
incoherent nuclear spin fluctuations. The question as to 
what really happens when an electronic moment tunnels, 
while it is strongly coupled to its nuclear spin, has not yet 
a clear answer. Contrary to the 3d group, hyperfine inter- 
actions are very large in 4/ elements. Diluted rare-earth 
ions in a nonmagnetic insulating single crystal are there- 
fore very suitable to study the possible entanglement of 
nuclear and electronic moments, when tunneling occurs. 
Our choice was the weakly doped rare-earth fluoride se- 
ries LiYi_ x R x F4, in which high quality single crystals are 
mainly investigated for applications in high-power laser 
diodes [O . Note that EPR spin-echo of magnetic tunnel- 



ing states have already been observed in a 1% Dy-doped 
crystal JlJ] . At higher concentrations, these crystals were 
used for phase transition studies of dipolar ordered mag- 
nets [ [l3|| . Among them, the holmium doped fluoride is 
a random, dipolar coupled system with an Ising ground 
state doublet (g e g « 13.3 fL4j , see also and refer- 
ences therein) and a pure isotope 1 = 7/2 nuclear spin. 
The magnetic properties of the Ising ferromagnet LHI0F4 
and spin glass LiYo.833Hoo.i67F4 have been studied by ac 
susceptibility |16|. In particular, enhanced quantum fluc- 
tuations, leading to a cross-over to a quantum spin glass, 
were evidenced in an applied transverse field |l7| . 

In this letter, we investigate a single crystalline 0.2% 
holmium doped LiYF4 at subkelvin temperatures. The 
isolated magnetic moments are weakly coupled by dipolar 
interactions (/io^dip ~ few mT) so that this very diluted 
insulator should exhibit a nearly single ion quantum be- 
havior, in which we are interested in continuity with our 
studies on molecular magnets. The crystal has a tetrago- 
nal scheelite structure with a C^h space symmetry group 
(74i/a), and the point symmetry group at Ho 3+ sites is 
S4, almost equivalent to T>2d (for LiHoF4, unit cell pa- 
rameters are a = b = 5.175 A and c = 10.74 A [p[). Be- 
cause of a very strong spin-orbit coupling, each magnetic 
ion of 165 Ho is characterized by its J = 8 ground state 
manifold (gj = 5/4), split by crystal field effects. These 
last give rise to a large uniaxial magnetic anisotropy, i.e., 
a high energy barrier hindering the magnetic moment 
reversal. However, we will see that quantum fluctua- 
tions due to significant transverse anisotropy terms dras- 
tically reduce this barrier. This effect was very much 
weaker in Mni 2 -ac (« 10% barrier reduction ||). Using 
the I J, M > basis and T)2d symmetry, the approximate 
Hamiltonian including hyperfine interaction writes 
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The a j, 7,7, and O™ are Stevens' coefficients and 
equivalent operators |19|j . Exact diagonalization of the 
136-dimensional Hamiltonian (1) was performed, using a 
set of crystal field parameters obtained by high-resolution 
optical spectroscopy: B\ = 273.9 K, B% = -97.7 K, 
Bl = -6.5 K, B\ = -1289.1 K, B% = -631.6 K jn). 

was taken as a fitting parameter of the measured res- 
onances (see below). J mixing, Jahn- Teller effect and 
nuclear quadrupole interaction are assumed to be negli- 
gible. We first show the results with Aj = in Fig. |l|. 
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1. (a): Energy levels vs average value of J z , in zero 
applied field, showing an Ising ground state doublet and a first 
excited singlet at « 9.5 K above, (b): Low-energy part of the 
Zeeman diagram. The first excited state (« 25 K below the 
next excited T2 singlet) defines an energy barrier of ~ 9.5 K 
hindering the magnetic moment reversal. 

The eigenstates transform as one of the four irreducible 
representations Ti r 2,3,4 of the S4 point group. Significant 
transverse crystal field terms B\0\ and B§0\ mix free 
ion states, with AM = ±4, so that eigenvectors are lin- 
ear combinations either of | ± 7 >, | ± 3 >, | T 1 > and 
|=F 5 > for r 3 . 4 ; I ± 6 >, I ± 2 >, I T 2 > and I T 6 > for 
r 2 ; or |± 8 >, I ± 4 >, |0 >, I T 4 > and I T 8 > for T ± . 
In Fig. [j]a), the calculated low-lying states within the 5 Is 
multiplet show a Ising ground state doublet noted as 
Itpf' >, while the first excited state, a T2 singlet, stands 
at k 9.5 K above (direct measurements give a 10 ± 1 K 
p5|). Fig. 0(a) also shows how the expected large bar- 
rier ~ 10 2 K is shortcut by large tunneling gaps due 
to transverse crystal field terms (emphasized by shaded 
areas between singlets belonging to the same represen- 
tation). A strong electronic level repulsion in the low- 
lying excited T2 states is clearly shown in Fig. 0(b) . This 
defines the energy barrier the magnetic moment has to 
overcome in order to reverse its polarization. At very low 
temperatures, the system should be equivalent to a two- 
level system with an energy barrier of rj 10 K. Actually, 



this picture is strongly modified when intraionic dipo- 
lar interactions are taken into account (Aj ^ 0). They 
lead to a more complex diagram in the electronic ground 
state, showing several level crossings for resonant values 
H„ (—7 < n < 7) in Fig. |^. The transverse hyperfine 
contribution ^Aj(J + I- + J-I+) induces some avoided 
level crossings between \ij)^,I z i > and ,I Z 2 >, with 
AI = \I Z 2 — I z \\, only when A//2 is an odd integer so 
that the two electronic low- lying states l^ > are cou- 
pled through nondegenerated excited electronic levels. 
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FIG. 2. Splitting of the electronic ground state doublet by 
the hyperfine interaction (Aj/ks ~ 38.6 mK; see below). The 
level crossings occur for resonant values of the longitudinal 
field H n ( — 7 < n < 7). Some are avoided level crossings and 
hyperfine levels repulsion is then induced by the electronic 
level repulsion in the excited states shown in Fig. 0(b) . 

Magnetic measurements were made at 0.04 < T < 1 K 
and for /lo-ff < 2 T, with a micro-SQUID magnetome- 
ter pfl] allowing field sweep rates up to 1 T/s. The 
crystal is first saturated in a large positive field ap- 
plied along the c— axis fJ,oH sa ,t ~ 0.3 T, and then the 
field H z is swept between ±7? sat . At slow field sweep 
rates, an isothermal process occurs, leading to staircase- 
like hysteresis loops at T < 200 mK, as shown in Fig. [|. 
These well defined steps come from quantum relaxation 
at avoided level crossings. At T = 50 mK, the lowest en- 
ergy level is mainly populated and magnetization steps 
are observed for — 1 < n < 3. The ground and excited 
tunnel splittings at H z = 0, being due to very small 
perturbations such as internal transverse fields, should 
be very similar. Therefore, quantum tunneling in zero 
field is mainly associated with the dynamics of the low- 
est avoided level crossing (the first excited crossing is at 
AE/ks — .9cffAiB-Hi/&B ~ 205 mK, assuming <? e ff ~ 13.3 
The amplitude of the next step, the resonance 
n = 1 at ^o-Hi = 23 mT, is about ten times larger, 
suggesting a larger tunnel splitting A. Indeed, Fig. g 
shows that the hyperfine induced tunnel splitting of the 
third excited avoided level crossing is large enough to 
render the barrier transparent (A 25 mK). The relax- 
ation time is thus simply given by thermal activation 
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t = T exp(2AE/kBT), with a long to because spin- 
lattice relaxation time T\ can be hours at very low tem- 
peratures and/or as a result of internal fields fluctuations. 
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FIG. 3. Hysteresis loops for v = 0.55 mT/s and at different 
temperatures. Inset: derivative of the magnetization normal- 
ized to M s , ddm/dH, at T = 50 mK and for v = 0.55 mT/s. 

In the same trend, the measured magnetization step 
ratio AM (n = l)/AM(n = -1) w 25 at T = 50 mK is 
approximately equal to the Boltzmann ratio, which con- 
firms that thermally activated quantum relaxation occurs 
at H z y£ 0. The barrier, essentially transparent due to 
this large splitting, becomes again finite out of resonance. 
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FIG. 4. Hysteresis loops in a constant transverse field at 
T — 50 mK and for v — 0.55 mT/s. A transverse field en- 
hances the quantum fluctuations in zero longitudinal applied 
field leading to a larger magnetization step. Inset: details of 
the Zeeman diagram around zero field. Thermally activated 
tunneling shows two possible channels over the first and the 
third, more efficient, excited avoided level crossings. 

The quantum relaxation is strongly enhanced by a con- 
stant transverse field, as a result of an increase of the 
tunnel splittings (see Fig. |^). In zero longitudinal field, 
the small tunnel splittings rapidly increase and hysteresis 
vanishes. A saturation of the magnetization at M ~ is 



observed in transverse fields larger than 100 mT, when 
the barrier is nearly transparent, and the small "over- 
shot" with an oscillation in M may be due to spin- 
phonon transitions. As expected for a large tunnel split- 
ting, sensitivity to a small transverse field is very weak for 
the resonance n = — 1. The inset in Fig. || shows dm/dH 
at T = 50 mK. The width of the resonant transitions is 
about hqAH = 2 — 3 mT which is expected from dipolar 
broadening. Similarly to molecular magnets, quasi-static 
fields due to dipolar interactions lead to a distribution 
of internal fields whereas field fluctuations, essentially of 
F~ nuclear spins, give homogeneous level broadening. 

A hysteresis loop measured at T — 50 mK for a much 
faster field sweep rate (v = 0.3 T/s) is shown in Fig. ||(a). 
A succession of equally spaced large and weak magneti- 
zation steps occur at fields H n , with — 14 < 2n < 14. 
The larger ones, with integer n, are associated with sev- 
eral equally spaced level crossings and the smaller steps, 
with half integer n, fall just in between when the levels 
are equally spaced (see Fig. ||). dm/dH is used to de- 
termine the H n values plotted in Fig. |^(a) inset. From 
the slope, we accurately obtain ^qH u = n x 23 mT. The 
electronic ground doublet is thus split by hyperfine in- 
teraction in eight doublets over an energy range of about 
1.44 K. We deduce Aj/k B ~ 38.62 mK, to be compared 
to Aj/Jcb ~ 40.95 mK |Q. The observed hysteresis loops 
depend sensitively on sample thermalization, showing 
that the spin-phonon system is not at equilibrium with 
the cryostat, leading to a phonon bottleneck pfj| , pT| , p2[ . 
At a fast field sweep rate v = 0.3 T/s, the system enters 
such a regime at T « 1 K (moderate sample thermal- 
ization) showing hysteresis without any magnetization 
steps down to T « 600 mK. When the field is swept 
back and forth, a stationary regime occurs and hystere- 
sis loops become nearly temperature-independent below 
a temperature T c (v) depending on sample thermalization 
(T c w 200 mK for v = 0.3 T/s). Below T w 600 mK, 
a nearly adiabatic process occurs, due to a much longer 
spin-lattice relaxation time Tj . The spin system becomes 
more and more isolated from the phonon bath, and en- 
ergy exchange between electronic and nuclear spins is 
only possible at fields H n . Equilibrium within the spin 
system is due to either quantum fluctuations at avoided 
level crossings (integer n) or to spin-phonon transitions 
and/or cross-spin relaxation, allowed by weak dipolar in- 
teractions, when energy levels are almost equally spaced 
(integer and half integer n) |^,^3|. Spin-spin interac- 
tions allow two additional steps for n = 8 and n = 9, 
at fields with equally spaced levels but no level cross- 
ing [Fig. ||(b) inset]. A small transverse applied field 
only increases the zero-field magnetization step, show- 
ing the weak effect of enhanced quantum fluctuations on 
hysteresis loops in this regime. Other resonances and 
small magnetization steps, dominated by cross-spin re- 
laxation, are not affected by a small transverse field, if 
small enough (fioHr ^ 20 mT). If the field sweep is sud- 



3 



denly stopped, the spin-phonon system exchanges energy 
with the cryostat and the magnetization relaxes toward 
the equilibrium curve. 
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FIG. 5. (a): Hysteresis loops at T = 50 mK and for 
v = 0.3 T/s. Several magnetization steps are observed for 
resonant values of the applied field fioHn ~ n x 23 mT [see 
inset; H n values are deduced from Fig. |(b)] . (b): Deriva- 
tive of the loop shown in (a) for a decreasing field. The two 
additional measured steps shown in the inset, for n = 8 and 
7i = 9, are associated with cross-spin relaxation only. 

The asymmetry of the envelope (not drawn) of the 
peaks in Fig. ||(b), showing that spins reverse mostly 
after field inversion, and the absence of constriction in 
the hysteresis loop near H z = 0, confirm the existence of 
small barriers, mainly in zero field (tunneling). 

In conclusion, we have shown that the quantum 
rotation of weakly coupled magnetic moments in 
LiYo.998Hoo.oo2F4 can be driven and monitored by hy- 
perfine couplings. At very low temperatures, when the 
field is slowly swept from negative to positive values, 
the coupled electronic and nuclear moments tunnels from 
\4>XiIzi > to \ip^,I Z 2 >• In a constant transverse field, 
the magnetization step, associated with incoherent tun- 
neling at the avoided level crossing in zero field, increases 
very rapidly. It saturates when the barrier is completely 
transparent. Details of hysteresis loops are in excel- 
lent agreement with the level structure of the electronic 



ground state doublet split by hyperfine interaction in six- 
teen states. At faster field sweep rates, additional magne- 
tization steps are observed and attributed to cross-spin 
relaxation and/or spin-phonon transitions in a phonon 
bottleneck regime. Very diluted Holmium doped LiYF4 
is thus a model system to study tunneling of an electronic 
moment strongly coupled to its nuclear spin. 
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